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We present experimental and computational studies of the propagation of internal waves 
in a stratified fluid with an exponential density profile that models the deep ocean. The 
buoyancy frequency profile N{z) (proportional to the square root of the density gradient) 
smoothly varies by nearly an order of magnitude over the fluid depth, rather than being 
■ constant or piecewise-constant as in prior studies. In addition to being nonuniform, the 

stratification is characterized by a turning depth Zc, where N{zc) is equal to the wave 
frequency uo and N{z < Zc) < uj. Internal waves reflect from the turning depth and 
become evanescent below. The energy flux below the turning depth is shown to decay 
exponentially with a decay constant given by /Cc, which is the horizontal wavenumber at 
the turning depth. The vertical velocity fields of the incoming and reflected waves above 
the turning depth agree within a few percent with a previously untested theory for a 
[ fluid of arbitrary stratification [Kistovich and Chashechkin, J. App. Mech. Tech. Phys. 

_C1< 39, 729-737 (1998)]. 
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1. Introduction 

In a stratified fluid, an internal wave transports momentum and energy as the prop- 
agating disturbance is restored by buoyancy forces. In the absence of rotation, internal 
waves satisfy the following dispersion relation 

sin(^(^)), (1.1) 



C/2 



Niz) 

where uj is the angular frequency of the wave, is the angle of propagation relative to the 
horizontal, N{z) = ^/— (g/po) (dp/dz) is the buoyancy frequency, g is the gravitational 
acceleration, po is a reference density and p{z) is the density profile. Approximately 
half of the i nternal wave energy in the ocean is produced by tidal flow over bottom 
topography (|Munk Wunschlll998l : IWunsch fc Ferraril l2QQ4V resulting in waves at the 



lunar semidiurnal tidal frequency ujm2 = 1.0452 x 10""^ rad/s, called internal tides. 

Oceanic buoyancy frequencies vary greatly, decreasing from the values in the thermo- 
cline t o become two orders o f magnitude smaller in the abyss where the water is well 
mixed. IPingree fc Newl ( 199ll ) found that deep in the Bay of Biscay the value of N{z) 



became so small that the internal waves became n early vertical (Eg . II. ip at the point 
of reflection from the bottom topography. Recently 'King et~al\ (|2Q12[ ) analyzed data for 
pressure, temperature, and salinity for thousands of locations throughout the oceans and 
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found that in many locations in the deep oceans the buoyancy frequency becomes even 
smaller than ujm2- Internal waves approaching a turning depth Zc where N{zc) = ujm2 
become vertical by (|l.ip and then reflect from the turning surface; below the turning 
depth the waves are evanescent (exponentially damped). 



Re c ent studies have examined i ntern a l waves in piecewise-linear stra tifications ([Sutherland fc Yewchuk 
20041 : iMathur fc Peacockl l2009l I2OIOI : iGregorv fc Sutherland! l2010[ ). Building on that 



work, we examine internal wave propagation in stratified fiuids with co ntinuously varying 

7V(z), as in the ocean. This physically important case was considered bv lKistovich fc Chashechkin 
(| 19981 ). who developed a linear theory that describes the propagation of internal waves 
in arbitrary stratifications, including those that contain a turning depth. 

We have conducted laboratory experiments and direct numerical simulations for a 
stratified fiuid with N{z) designed to mimic ocean conditions, including a turning depth. 
Our laboratory and numerical results are found to agree well, and our results for refiected 
and evanescent waves agree within a few percent with the theoretical predictions of 
Kistovich and Chashechkin. 

The re st of the paper is organize d as follows. In ^we summarize the results of the 
theory of iKistovich fc Chashechkinl ^1998) for wave refiection from a turning depth. In 
^ we describe our numerical and experimental methods. In ^ we present our results, 
the consequences of which are discussed in §5] 



2. Kistovich-Chashechkin theory 

The theory of IKistovich fc Chashechkinl (19981) describes the propagation of internal 
wave beams in arbitrary stratifications. The theory accounts for the role of viscous dissi- 
pation, diffusion and the interaction of an internal wave beam with a turning depth. The 
density profile is assumed to be po{z) = poo [1 + so{z)]i where so is the reduced salinity 
and poo is a reference density. The internal wavefield is assumed to be two-dimensional 
and to have a time-dependence of the form exp(— icjt) (which is henceforth omitted). 
Then the linear Navier- Stokes equations in the Boussinesq approximation become 

1 dp . . 1 dp 

— luju = — — h h'Au^ —luw = — 

Poo ox Poo oz 

dso ^ . du dw 
-lujs ^ w—— = DAs, + ^ 

az ox oz 

where u is the horizontal velocity, w the vertical velocity, p and s the variable pressure 
and salinity, u the kinematic viscosity, D the salt diffusivity, g gravitational acceleration 
(in the direction — i), and A = d'^/dx'^ + jdz^ is the two-dimensional Laplacian. 

The internal wavefield may be described by the vertical displacement /i(x, z\ which is 
related to the vertical velocity by w = —\uoh. Solutions may be sought in the form 



0, 



(2.1) 
(2.2) 
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where k is the wave number. 

Asymptotic solutions for / may be obtained if N{z) varies slowly compared to the 
length scales describing the internal wave beam. Given their assumptions, Kistovich and 
Chashechkin obtain 
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where A{k) describes the spectral properties of the wave source located at (xo^zo), v 
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V + D, = 1 (—1) corresponds to beams propagating down (up) and to the right, and 

iV2(2) -w2 



(2.5) 



Thus, given A{k) and N{z)^ the theory of Kistovich and Chashechkin predicts the vertical 
velocity field everywhere in space. 

In addition to examining the propagation of internal wave beams in arbitrary strati- 
fications, Kistovich and Chashechkin derived the following expressions for the refiection 
of an internal wave beam from a turning depth Zc, defined as N{zc) = uo and N{z) < uo 
for all z < Zc'. 
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where fi describes the incoming wave beam for x < Xc and 2: > Zc, /r is the refiected 
wave for x > Xc and z > Zc and ft is the transmitted (evanescent) wave for z < Zc. Given 
the spectral amplitudes A{k), we have 



B{k) = A{k)e-''^'e 



i0Opi7r/2 



C{k) = A{k)e-''^'e 
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where 



k = k f 



7- 



k^ (1 + 7')' 



26^7 



dz\ 



(2.9) 
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which takes into account the phase winding and viscous decay as the beam propagates 
from the source at (xq^zq) to the turning depth at (xc^Zc). Therefore, only N{z) and 
A{k) are required to fully characterize the vertical velocity field. 



3. Methods 

3.1. Direct Numerical Simulations 

We numerically simulate internal wave propagation in nonli near stratification s using 
CDP 2.4, which is a large eddy simulation code developed bv iHam fc laccarinol ([2004). 
By disabling all subgrid-scale modeling, we perform direct numerical simulations of the 
incompressible Navier-Stokes equation. The code is a parallel, unstructured finite- volume- 
based solver that is modeled after the algorithm presented by Mahesh et al. ([2004) . 



Second-order accuracy in space and time is achieved by using a fractiona l step time 
marching scheme and multiple implicit schemes for the spatial operators (' Ham et~al\ 
2006h . We perform two-dimensional simulations in the Boussinesq approximation. These 



approximations are justified by the small density variation and predominantly two- 
dimensional fiows examined in the experiments. CDP solves the following equations for 
the pressure p, density p and velocity v = {u^ w) in the (x, z) directions: 

dv/dt + {v ' V)v = -Vp/po - gp/poz + i^V^v + F/po, (3.1) 
V'V = 0, dp/dt^{v'V)p = DW'^p, (3.2) 
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where po = 1000 kg/m^ is a reference density, g is the gravitational acceleration, v is the 
kinematic viscosity, F drives the internal waves, as described below, and D = 10~^^ m^/s 
is the density diffusion coefficient. The value of D is so small that salt diffusion in the 
simulations is negligible; this is an accurate representation of the experiments. 

Internal waves are generated by th e forcing term . F show n in Eq. (13.11), which is simila r 



to that used bv lSlinn fc RilevI (|l998[ ). Ijavam et all (|l999h and iRodenborn et all (|201lh : 



F{x, = (V X V X ^z) sm{ujt), (3.3) 
^{x,z) = ^oexp - (^a{x - xq)'^ ^2b{x - xo) {z - zo) ^ c{z - zo)'^^ , (3.4) 

where a = cos'^{f3)/2al + sin^(/3)/2a^, b = sin(2/3)/4cr^ - sin(2/3)/4cr^, c = sin^(/3)/2cr^ + 
cos^(/3)/2<j^, uj is the internal wave frequency, the maximum forcing amplitude, 
(xq^zq) the spatial center, P the counterclockwise rotation of the elliptical Gaussian 
forcing profile, and ax {(Tz) is the horizontal (vertical) standard deviation. The angle /3 is 
chosen so that the forcing profile is aligned with the local internal wave propagation angle 
0. The oscillation frequency is cuq = 0.628 rad/s for all of the simulations and experiments 
except for Run 12 where uj = ujq/A (see table [1]). The time step is At = 7r/4000cc;, which 
yields 8000 time steps per period. We run the simulations for 40 periods to ensure that 
the system has reached a steady-state. The forcing amplitude and standard deviations 
ax and are varied to modify the properties of the wave beam, as described in table [H 

Two computational grids are generated using Pointwise Gridgen. Grid I is designed for 
direct comparison of the numerical simulations with the experiments. The domain has 
the same dimensions as the experimental tank, namely < x < 0.9 m and < z < 0.6 m. 
The numerical wavemaker is centered at (xq = 0.15 m, zq = 0.50 m) in both the simula- 
tions and experiments. While the grid is structured everywhere (composed of rectangular 
control volumes), the horizontal and vertical resolutions vary smoothly throughout the 
domain to take into account the variable shear as the beam propagates in the nonlinear 
stratification. Near the wavemaker the wave beam is nearly horizontal, resulting in strong 
vertical shear, while at the turning depth the beam propagates vertically, necessitating 
increased horizontal resolution. At the wavemaker center we use a horizontal resolution 
of Ax = 2.5 mm, while the vertical resolution is Az = 0.4 mm. The resolutions near the 
reflection from the turning depth at (xc, Zc) are Ax = 0.4 mm and Az = 1.5 mm. Grid I 
is composed of approximately 5.5 x 10^ control volumes. No-slip boundary conditions 
are enforced at x = 0, z = 0, and x = 0.90 m. The top surface at 2: = 0.60 m is a 
free-slip, impenetrable boundary chosen to emulate the free surface in the experiments. 
The internal waves are damped for x > 0.8 m by adding a drag force equal to — 0.002 v 
to Eq.( I3.ip to mimic the absorbent material present in the experiments. 

Grid II is designed to compare th e numerical simulations with the theoretical pre- 
dictions of iKistovich fc ChashechkinI ()l998[ ), which are described in ^ The grid spans 



< X < 6 m and < 2; < 2 m. The wavemaker is centered at {xq = 1.0 m, = 1.5 m). 
As with Grid I, the spatial resolution varies smoothly throughout the domain to account 
for the variation in the beam angle as it propagates. The horizontal resolution in the 
vicinity of the wavemaker is Ax = 15 mm while the vertical resolution is Az = 2 mm. 
Near the reflection at the turning depth the horizontal resolution is Ax = 4 mm and the 
vertical resolution is Az = 6 mm. Grid II has a total of ~ 2.5 x 10^ control volumes. 
No-slip boundary conditions are applied at all of the boundaries. To reduce finite-size 
effects, a drag force of the form —0.02v is added to Eq. 13. II for x < 0.6 m, x > 5.5 m and 
z > 1.75 m to damp the wave beams before they reach the boundaries. 

The density profiles are chosen to be scaled models of the stratification in the deep 
ocean. In particular, we are interested in examining the dynamics of internal waves at 
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Figure 1. (a) The experimental density profile (black squares) agrees well in the domain of 
interest with the density profile given by Eq. (|3.5p and used in the simulations (red line), (b) 
The experimental wave maker consists of five acrylic plates in a box (right) that are driven by 
a rotating camshaft (left). Each successive cam is rotated by 45°, which results in a half-sine 
wave configuration. 



turning depths analogous to those recently discovered by iKing et al\ (l2Q12h . As such, 
all of our density profiles are exponentials with a turning depth at Zc, where N{z) < u 
for z < Zc. However, the results presented here do not depend upon the stratification 
being exponential in form; the crucial aspect is the existence of a turning depth. The 
stratification used for the experimental domain covered by Grid I is shown as the red 
line in Fig. [H and is given by 

3 



Pexpt = 1061 - 10exp(3.12z) kg/m'', ^ z ^ 0.6 m. 



(3.5) 



The depth-dependent buoyancy frequency A/expt(^) has a turning depth at Zc = 0.082 m. 
The density profile used for all simulations with Grid II, except Run 12, is given by 



1244- 5.75 exp(2.5z), 
1916 -611z, 



^ z ^ 1.5 m 
1.5 < z ^ 2 m. 



(3.6) 



The resulting buoyancy frequency profile varies by a factor of 6.5 from top to bottom and 
has a turning depth at Zc = 0.424 m. The density profile for Run 12 has N{z) reduced 
by a factor of 4 for all depths. 

3.2. Experimental Methods 

Experiments are performed in a glass tank with < x < 0.90 m, < y < 0.45 m and 
< z < 0.60 m. The exponential gradient given b v Eq . (|331) is generated by utilizing the 
generalize d version of the double-bucket method (jFortuinlll960l : lo"sterll 19651 ) described by 



Hilll (2002). The density profile is measured by extracting fluid samples at known heights, 
which are then measured with an Anton Paar density meter. An example experimental 
stratification is shown as the black squares in Fig.dja). The fluid at the top {z = 0.59 m) 
has p = 1000 kg/m^. The experimental density profile agrees well with the stratification 
used in the simulations (red line and Eq. (j3.5p ) for 0.1 < z < 0.5 m. 

Internal wave beams are generated using a wav e maker similar to the one developed by 
3stiaux et all (120071 ). Our wave maker (Rodenb orn et a/.ll201l[ ) consists of five acrylic 
plates with dimensions 150 mm x 150 mm x 6 mm that are separated by 2.5 mm and 
placed within a parallelpiped open-sided box, as shown in Fig. [TJb). A helicoidal rotat- 
ing camshaft oscillates the plates periodically with the direction of rotation setting the 
internal wave phase velocity vector. The eccentricity of the camshaft (12.7 mm) deter- 
mines the oscillation amplitude. The generated wave beams are nearly two-dimensional, 
particularly along the center of the tank {y = 0.225 m) where all of our measurements 
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Figure 2. Snapshot of the vorticity (color) for a wave beam generated in the simulations at 
(x = 1 m, z = 1.5 m) after 40 periods of oscillation. The stratification (Eq. ()3.6|) ) has a turning 
depth at Zc — 0.424 m (dashed line). 



are taken. The wave maker is centered at (xq =0.15 m, = 0.50 m) and oriented such 
that the plates are parallel to the the wave beam angle at z = zq. The internal wave 
motions are dissipated by a fi ber mesh for x > 0.88 m to reduce finite-size effects. 

Particle image velocimetry (lAdrianlll99l"h is used to obtain the two-dimensional velocity 
field V = (ii, w) in the vertical plane given hy y = 0.225 m. Hollow glass spheres with 
diameters S < d < 12 /im and densities in the range 1050 < p < 1150 kg/m^ are used 
as seed particles. The tracers are illuminated by a 5 mm thick laser sheet produced by 
a 532 nm wavelength laser (2.5 W). Two 12-bit CCD cameras with 1296 x 966 pixel 
resolution capture orthogonally scattered light in a 194 x 145 mm^ region, 40 times 
per period {At = 0.25 s). The instantaneous v elocity fields are determined using the 
CIV algorithm developed bv lFincham fc Delercd (j2000) and are interpolated to a regular 
50 X 50 grid with spatial resolutions of Ax = 3.9 mm and Az = 2.9 mm. 



4. Results 

An example snapshot of a wave beam propagating in the stratification given by Eq. 
(|3.6p is shown in Fig. [2] (see movies in online supplementary material), with vorticity 
Q = V X V indicated by color and the turning depth z = Zc given by the dashed hue. As 
can be seen, the beam generated at (x = 1.0 m, z = 1.5 m) propagates down and to the 
right before being reflected at the turning depth, causing the energy to move up and to 
the right. Only a small fraction of the energy flux penetrates below the turning depth. 

We compare the velocity and vorticity fields of the numerical simulations with the 
experimental measurements to verify the computations, as illustrated by the incoming 
and reflected beam profiles in Fig. [3l Horizontal profiles of the vertical velocity are well 
described by modulated Gaussians of the form 

w{x,z,t) = Ao{z) COS {ko{z)x - ujt)ex.p I- [x - p{z)] /2a{z) j [j{zo)/j{z)] ^ , (4.1) 

where Aq is the velocity amplitude, ko the wavenumber, ja the beam center, a the beam 

1/2 

width, and the term [y{zo)/'y{z)] ' accounts for the increase (decrease) in the vertical 
velocity as the beam propagates into weaker (stronger) stratification. Therefore, even 
though the vertical velocity may increase by geometrical effects as a beam propagates 
into weaker stratifications, the velocity amplitude monotonically decreases owing to 
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Figure 3. Snapshots of the vertical velocity (black) and vorticity (red) horizontal cross-sections 
centered at (a) (/x = 0.30 m, z = 0.40 m) for the incoming wave beam and (b) (/x = 0.76 m, 
^ = 0.37 m) for the reflected wave beam. Laboratory measurements are given by the data points 
and results from the simulations are shown as solid lines. 



viscous dissipation. We find that the experimental and computational vertical velocity 
and vorticity profiles are in remarkable agreement, as Fig. [3] illustrates. 

4.1. Evanescent internal waves 

Below the turning depth Zc, internal waves be come evanescent with an intensity that 
decays exponentially with a decay constant a ( Gostiaux et fl/.|[2006[ ). We compute the 
intensity of the evanescent wave beams by horizontally integrating the vertical flux at 
different depths, 



$(z) 



(4.2) 



where and are the fluctuations in the pressure and vertical velocity, respectively. 
Figure Hl^a) shows example results for the vertical flux below the turning depth. The 
distance from the turning depth Zc — z is scaled by the horizontal wavenumber at the 
turning depth kc = A^o(^c), which is given near each curve. Below the turning depth, the 
flux decays exponentially, except for large Zc — z where bottom-boundary interactions 
arrest the decay. 

Our measurements of the decay constants a agree well with the expected values a = kc 
(see Fig. IH^b) and the summary in table [2j). Therefore, in analogy with evanescence in 
other contexts, we conclude that the flux is well described by 

-kc{zc — z) 



for z < Zc 



(4.3) 



4.2. Internal wave reflection: comparison with theory 

The theory of 'Kist ovich fc Chashechkinl ( 1998 1) described in ^provides predictions with 
no free parameters for the behavior of the vertical velocity fleld of internal waves prop- 
agating in nonlinear stratiflcations, even those with a turning depth. To test these pre- 
dictions, we compare horizontal cross-sections of the vertical velocity fleld computed by 
our numerical simulations to the theoretical predictions. Beam proflles are flt to the 
form (j4.ip to determine the amplitude Ao{z)^ wavenumber ko{z)^ beam center jj.{z) and 
beam width cr(z), as the beam propagates towards the turning depth, reflects, and then 
propagates away. 
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Figure 4. (a) Results from the numerical simulations for the horizontally integrated flux 
defined by Eq. ()3.4|) . for evanescent waves beneath the turning depth Zc. The horizon- 
tal wavenumber at the turning depth, /cc, is indicated next to the corresponding curve, (b) The 
decay constant a from simulations (black circles) is compared to the horizontal wavenumber at 
the turning depth, kc (red line). 



To use the theoretical expressions Eqs. (j2.6p and (|2.7p for the incoming and reflected 
wave beams, we must know the stratification N{z) and the spectral amplitudes of the 
wave source A{k). The buoyancy frequency profile may be determined using Eq. ()3.6p . 
We determine A(k) by fitting a horizontal cross-section of the vertical velocity near the 
wave source to Eq. ()4.ip . which we then Fourier transform analytically to determine 
the spectral amplitudes. We then compute the predicted vertical velocity field for the 
incoming and reflected wave beams using Eqs. ()2.6p and (j2.7p . 

The numerical simulations and theoretical predictions are compared by measuring 
the velocity amplitude Aq{z)^ using the same values of {A:o(^), /i(^), and (j{z)} for the 
simulations and theory. The simulation results for Aq^z) as a function of the geometric 
length from the turning depth (L — Lc) are compared to theory in Fig. [5l where 

L{z)= r ^/l + j^{z')dz', (4.4) 

J Zo 

is the distance from the wave source measured along the centerline of the beam, and 
Lc = L{zc). The four curves in Fig. [5fa) correspond to different viscosity values, as 
indicated. Figure Efb) compares cases with the same viscosity but different values of the 
source wavenumber, source amplitude, or wave frequency, as indicated. In all cases that 
we have studied, the fully-constrained theoretical predictions agree (on average) within 
1.5% of the numerical simulations, as summarized in table [2j 



5. Discussion 

We have presented laboratory experiments and direct numerical simulations of inter- 
nal wave beams propagating in a fluid where the buoyancy frequency N{z) decreases 
exponentially as the bottom boundary is approached, as is often found in the abyssal 
ocean. We have examined wave propagation in a fluid wi th a turning dep th Zc, below 
which N{z) < uj^ diS occurs at many locations in the ocean ( King et al. 2Q12[ ). The exper- 
iments and simulations agree within a few percent with the predictions of the previously 
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Figure 5. The vertical velocity amplitude ^0(2;) in the simulations (data points) is compared 
with the fully constrained theoretical predictions of Kistovich and Chashechkin (solid lines) (a) 
for varying viscosity values, and (b) for fixed viscosity (5 x 10~^m^/s) and other parameters 
given by {u, /co(^o), (t{zo)): Run 6, (0.63 rad/s, 18.7 m"\ 0.127 m); Run 9, (0.63 rad/s, 22.3 
m-\ 0.144 m); Run 12 (0.16 rad/s, 22.3 m-\ 0.144 m). 



untested linear theory of iKistovich fc Chashechkin ( 1998 ) for internal wave beam prop- 
agation in a fluid with arbitrary stratification. Thus if the stratification N{z) and the 
internal wave beam spectral amplitude A{k) are known, the velocity field of a propagat- 
ing internal wave beam can be predicted with no free parameters, as long as topography 
and nonlinear interactions are negligible. 

The horizontally-integrated flux below a turning depth decays exponentially with a 
decay constant a = kc-, where kc is the horizontal wavenumber at the turning depth. This 
decay may greatly limit the role of internal tides beneath turning depths. For example, a 
wave beam with a wavelength A ~ 300 m incident upon a turning depth will produce an 
evanescent beam with a horizontally-integ rated flux that dec ays by an order of magnitude 
every ~ 100 m beneath the turning depth. iKing et al\ ( 2012[ ) found that in some locations 
turning depths occur at depths greater than 1 km above the ocean bottom. 

Another consequence of the existence of locations with turning depths well above the 
ocean floor is the effect on internal wave generation by tidal flow over topography in 
a region where internal waves are evanescent. If the turning depth is near the ocean 
floor, the internal waves could tunnel through the fluid beneath the turning depth, but 
in general internal tides generated below a turning depth will be greatly reduced in 
intensity by the exponential damping. Thus future studies should examine the effect of 
turning depths on internal tide generation for a wide range of parameters (e.g., turning 
depth height, N{z)^ tidal frequency). 
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Run 


cj (rad/s) ly 


(xlO-^ mVs) 


Ao{zo) (mm/s) 


ko{zo) (m ^) 


a{zo) (m; 


1 


0.63 


1 


0.29 


24.1 


0.109 


2 


0.63 


2 


0.27 


23.8 


0.111 


3 


0.63 


5 


0.22 


22.0 


0.124 


4 


0.63 


10 


0.18 


20.7 


0.135 


5 


0.63 


1 


1.37 


23.2 


0.123 


6 


0.63 


5 


1.13 


22.3 


0.127 


7 


0.63 


5 


3.38 


18.4 


0.157 


8 


0.63 


5 


0.12 


14.6 


0.176 


9 


0.63 


5 


0.16 


18.7 


0.144 


10 


0.63 


5 


0.25 


25.5 


0.117 


11 


0.63 


5 


0.09 


25.9 


0.137 


12 


0.16 


5 


0.33 


16.5 


0.189 



Table 1. Parameters describing the direct numeric al simulations perform e d wit h Grid 
II and psim, used to compare to the theory of iKistovich ChashechkinI (|l998l ). Hor- 
izontal cross-sections of the vertical velocity field near the wave source are given by 
w{x,t) = Ao(zo) cos {ko{zo)x — uot) e:><.p [—{x — /j{zo))'^ /2a (zo)], where Aq is the velocity am- 
plitude, cj the angular frequency, ko the wavenumber and a the beam width. The kinematic 
viscosity u is constant throughout the simulation domain. 



Run ^{zc) 


(10-^ kg/s^) 


kc (m-^) 


a (m-^) 


S{%) 


1 


-3.82 


22.2 


22.7 


0.48 


2 


-4.15 


21.0 


22.1 


1.08 


3 


-3.63 


18.6 


20.3 


0.96 


4 


-2.52 


17.0 


18.2 


2.08 


5 


-107 


22.2 


21.8 


2.28 


6 


-92.6 


18.6 


20.2 


0.50 


7 


-1721 


18.2 


17.9 


0.74 


8 


-1.81 


14.6 


15.5 


3.42 


9 


-2.70 


17.0 


17.7 


0.86 


10 


-3.89 


20.9 


20.0 


2.77 


11 


-0.91 


23.7 


25.2 


2.20 


12 


-4.90 


14.0 


15.3 


1.39 



Table 2. Summary of results presented in Figs. [H and [5l where ^{zc) is the horizontally-inte- 
grated vertical flux at the turning depth, kc is the horizontal wavenumber at the turning depth, 
a is the decay constant describing the exponential decay of the vertical flux below the turning 
depth, and 6 is the mean abs olute percent difference b etween our measurements of Ao(z) and 
the theoretical predictions of IKistovich ChashechkinI (jT998|). 



